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Abstract 



Recently it has become clear, that a key issue in quantum computation is understanding how interaction with 
the environment, or "decoherence" , effects the computational power of quantum computers. We adopt the standard 
physical method of describing systems which are interwound with their environment by "density matrices", and 
1 ' within this framework define a model of decoherence in quantum computation. 

Our results show that the computational power of decohered quantum computers depends strongly on the amount 
of parallelism in the computation. We first present a simulation of decohered sequential quantum computers, on 
a classical probabilistic Turing machine, and prove that the expected slowdown of this simulation is polynomial in 
' time and space of the quantum computation, for any non zero decoherence rate. Similar results hold for Quantum 
0^ ■ computers that are allowed to operate on logarithmic number of qubits at a time. 

For decohered quantum circuits (with local gates), the situation is more subtle and depends on the decoherence 
rate, r\. We find that our simulation is efficient for circuits with decoherence rate i] higher than some constant 
rji, but exponential for a general (random) circuit subjected to decoherence rate lower than some constant 772. The 
\Q , transition from exponential cost to polynomial cost happens in a short range of decoherence rates. We use computer 
Q\ experiments to exhibit the phase transitions in various quantum circuits. 

Oh! 

h— > ■ 1 Introduction 

Quantum Turing Machines^, and Quantum Circuits j?], [H| challenge the so called "polynomial Church thesis" 
H ' which asserts that "randomized Turing machines can simulate with polynomial slowdown any computational de- 
vice." In particular Shor's quantum factoring algorithm p6| provides within this theoretical framework, an efficient 
solution, to a problem for which no classical polynomial algorithm is known. 

It is yet unclear whether and how quantum computers will be physically realizable,^, |j| but as any physical 
■ system, they in principle will not be ideally isolated, and will to some extant interact with their environment. Such 
interaction causes the state of the computer to be interwound, or entangled, with the state of the environment, 
a process called "decoherence" [[l8| [I?], Q . A real challenge to the "polynomial church thesis" would come from 
decohered quantum computers. Early works Q showed that the effects of decoherence on the quantum computation 
can not be ignored, and that decoherence may limit the applicability of quantum algorithms. Thus, understanding 
the computational power of decohered quantum computers is essential. 

Decoherence takes us out of the standard model of quantum computers described by pure states, as the state of 
a decohered computer is in general a mixed state. Mixed states are used when there is only partial knowledge about 
the system, due to the fact that the system is interwound with it's environment. Such a state can be represented 
as a probability pk for the system to be in the pure state \ctk). This description is not unique and instead we use a 
unique representation of mixed states called density matrices fll|| , to describe the state of the quantum computer^). 

In order to incorporate decoherence into the model, we need to add a clock to the circuit since the process of 
errors is dynamic and depends on time. A circuit is associated a timing which indicates the time step for each 
gate to operate. To add decoherence, the following assumptions are made: (1) Each qubit decoheres independently 
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(single qubit faults), and (2) No decoherence takes place inside the gates. Each qubit decoheres, or undergoes a 
fault with probability rj per step, and r\ is referred to as the decoherence rate. The triplet of a quantum circuit, a 
timing of the circuit, and a decoherence rate construct together what we call a quantum medium, which operates 
subjected to different types of single qubit faults. The list of times and places where faults had occured, namely 
the fault path, is random, and naturally, the function that the medium computes is the weighted average over 
the outputs for each possible fault path. This model for decoherence is richer than it seems at first glance: it is 
equivalent to noise caused by weak interactions with the environment, and, to a model in which there are many 
possible kinds of faults that can occur. Our results therefore apply to these kinds of decoherences too. 

In this paper we exhibit upper bounds on the computational power of quantum mediums subjected to faults. As 
a first observation, note that regarding upper bounds, it is enough to restrict the discussion to faults of a specific 
type, which we choose to be collapses of a single qubit, or in other words: measurements of a single qubit in some 
basis. It suffices to deal with collapses because of two reasons: The first is that if a quantum medium claims 
to compute a function fault tolerantly, it is supposed to compute it when subjected to any fault. The second, 
and maybe more significant reason, is that it is reasonable to assume that every physical realization of quantum 
computers will be subjected to some extent to collapses . 

We present a simulation of a quantum medium subjected to single qubit collapses, on a classical probabilistic 
Turing machine. The simulation chooses, with the appropriate probability, a fault path. Then it keeps track of 
the development in time of the density matrix of the medium according to this path, and gets a final density 
matrix. The simulation imitates measurements on this density matrix and outputs the result. Since the algorithm 
chooses the path with the correct weight, the distribution over the results will be as that of the medium. How 
does the simulation keep track on the development of the density matrix, which is exponential in size? Trying to 
be efficient, the simulation presents this matrix as a tensor product of smaller density matrices. The input string 
is described as a tensor product of the density matrices of each qubit, but if no faults occured, the computer can 
soon entangle all the qubits together, with no way to view their exponential density matrix as a tensor product of 
smaller matrices. Collapses prevent this from happening, since after a qubit collapses to one of it's basic states it 
is no longer entangled to the other qubits, and can be described by it's own density matrix. 

We are interested in the expected cost of the simulation. The results emphasize the significance of the amount 
of parallelism in the computation. To gain some intuition, one can view the decohered quantum computation as a 
struggle between two forces: The computation, which tries to use the exponential dimension of the Hilbert space, 
by entangling as many qubits as it can, and the decoherence, which destroys this entanglement by collapses. ^From 
this point of view, a sequential computer obviously can not win the battle - each time step it applies only one 
local operation, but 9{n) qubits collapse. The picture in mediums with general parallelism is different, since both 
forces that compete have linear power: 6{n) gates can be applied each time step, and the decoherence collapses a 
percentage of the n qubits. Since the competition here is even, the question which force wins in this struggle is 
delicate and requires careful consideration. 

Our results show that: 

• The expected slowdown of our simulation applied on mediums with 0(log(n)) parallelism, with any 
non zero decoherence rate, is polynomial in the number of qubits and time steps. This holds also for multi head 
Quantum Turing machines that are allowed to operate on 0(log(n)) qubits each time step. 

The mathematical techniques involved in the case of mediums with general parallelism are more sophis- 
ticated. In our analysis we assume that the gates are of fan-in not bigger than two[[|, but the results can be 
extended to any constant fan-in: 

• High decoherence rate: The simulation of mediums with decoherence rate higher than some constant is 
efficient. Single-qubit-gates can not entangle different qubits, and therefore do not participate in the "struggle" 
against decoherence. Hence if two-qubits-gates are applied "rarely", the cost is polynomial already at lower 
decoherence rates. 

• Low decoherence rates: The cost of simulating a random circuit is exponential, for decoherence rates lower 
than some constant. 

• Phase transition in the cost: Results regarding phase transitions (To) imply that the cost of simulating a 
random circuit transforms from exponential to polynomial in a very short range of decoherence rates. As for other 
quantum circuits, if the cost transforms from exponential to polynomial this transition is sharp (by extension of 
JTc[ , private communication). 

We used computer experiments to show that in the case of a random circuit, the phase transition occurs at 
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?7o ~ 0.63. We also looked at a medium which is a one dimensional array of qubits, where gates are applied only 
on nearest neighbors, and found that in this case the transition occurs already at 77 ~ 0.50. 

Organization of paper: In section 2 we recall the basic definitions of quantum circuits and mixed states. 
In section 3 we define the model of quantum mediums with faults. Section 4 is devoted to the simple case of 
simulating sequential mediums. In section 5 we concentrate on simulating mediums with general parallelism. 
Section 6 discusses consequences of this work. 

2 Quantum circuits with mixed states 

In this section we recall some physical definitions of quantum systems of n two-state particles, Hilbert space, 
computational basis, pure states, mixed states and density matrices jn|. We then describe the model of quantum 
circuits with mixed states, defined in 

2.1 Pure states and mixed states 

Pure states: We deal with systems of n two-state quantum particles, or "qubits". The pure state of such a 
system is a unit vector, denoted \a), in the Hilbert spaced] C 2 " i.e. a 2™ dimensional complex space. We view 
C 2 as a tensor product n two dimensional spaces, each corresponding to a qubit: C 2 = C 2 <g> ... <£> C 2 . As a 
basis for C 2 , we use the 2™ orthogonal basic states: \i) = \ii) ® | ^2) - - - - ® |*n),0 < i < 2™, where i is in binary 
representation, and each L- gets or 1. Such a state corresponds to the j'th qubit being in the state A pure 

state I a) eC 2 is a superposition of the basic states: \a) = Y^i=i c i|*)i with J2i=i l c «| 2 = 1- l a ) corresponds to the 
vector v a = (cj.,C2, ...,C2»). uj,, the complex conjugate of v a , is denoted (a\. The inner product between \a) and 
|/?) is (a\/3) = (va,Vp). The matrix v^vp is denoted as An isolated system of n qubits develops in time by 

a unitary matrix]^] of size 2" x 2": \a(t 2 )) — U\a(ti)). A quantum system in C 2 can be observed by measuring the 
system. An important measurement is a basic measurement of a qubit q, of which the possible outcomes are 0, 1. 
For the state \a) — J2i=i C «N); the probability for outcome is po = J2i i\ =0 l c *l 2 ano - the state of the system will 
collapse to \(3) = ^ _ Ci\i), (the same for 1). In general, an observable O over C 2 is an hermitian^] matrix, of 

size 2™ x 2". To apply a measurement of O on a pure state \a) £ C 2 . write \a) uniquely as a superposition of unit 
eigenvectors of O: \a) = 2~^j Cj|oj), where \o{) have different eigenvalues. With probability \ci\ 2 the measurement's 
outcome will be the eigenvalue of |oj), and the state will collapse to \oi). A unitary operation U on k qubits can be 
applied on n qubits, n > k, by taking the extension U of U, i.e. the tensor product of U with an identity matrix 
on the other qubits. The same applies for an observable O to give O. 

Mixed states: A system which is not ideally isolated from it's environment is described by a mixed state. There 
are two equivalent description of mixed states: mixtures and density matrices. Two different density matrix can be 
distinguished by a measurement, and two systems with the same density matrix are quantumly indistinguishable. 
In contrast, the mixture description is not unique: different mixtures may yield the same density matrix. We use 
density matrices in this paper. 

mixtures: A system in the mixture {a} = {p^, \ctk)} is with probability pk in the pure state \ctk)- The rules of 
development in time and measurements for mixtures are obtained by applying classical probability to the rules 
for pure states. A unitary matrix U transforms a mixture {pk, \&k}} to {pk,U\ak}}- To apply an observable O 
on the mixture {a} = {pj, \cxj}}. write each pure state as a sum of eigenvectors of O with different eigenvalues: 
\aij) = J2i c j\°j)- The probability to get an eigenvalue A.; is Pr(Aj) = J2jPj\ c j\ 2 - ^ n the resulting mixture, 
under the condition the result was A^, each pure state has collapsed to it's part that is consistent with A^, and the 

conditioned probabilities are computed classically: {a}\i = \ ril — 1 1 })/- O ne can a ^ so a PPly an unconditioned 

measurement, i.e. not condition on the outcome of the measurement but instead average on the possible outcomes. 
This gives the mixed state O o {a} = {Pr(Xi), {a}^^ }. 

Density matrices: A density matrix p on C 2 is an hermitian positive semi definite complex matrix of dimentions 
2" x 2", with tr(p) = 1. A pure state \a) = Ci\i) is associated the density matrix pi a \ — \a)(a\ i.e. p\ a \(i,j) = 

1 A Hilbert space is a vector space with an inner product 

2 Unitary matrices preserve the norm of any vector and satisfy the condition U~ 1 = C/t 
3 An hermitian matrix H satisfies H = H* 
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CiC*. A mixture {a} — {pi, |a/)}, is associated the density matrix : p{ Q j = pip\ a ,) • The operations on a density 
matrix are defined such that the correspondence to mixtures is preserved. If a unitary matrix U transforms the 
mixture {a} = {pi, |a/)} to {/?} = {p u U\a{)}. Then p {/3} = Y^iPiU\ai)(ai\W = Up^W. Applying a measurement 
of an observable O on p, which is written in a basis of eigenvectors Vi of O. gives, the probability for an outcome A 
is the sum of the diagonal terms of p, which relate to the eigenvalue A: pr{\) — Pvi,vi$(K = A), conditioned 

that the outcome is the eigenvalue A, the resulting density matrix is 0\ o (p), which we get by hrst putting to zero 
all rows and columns in p, which relate to eigenvalues different from A, and then renormalizing this matrix to trace 
one. Without conditioning on the outcome the resulting density matrix will be O o (p) = Pr(\k)0\ k o (p). 
which differs from p, only in that the entries in p which connected between different eigenvalues are put to zero. 
Given a density matrix p of n qubits, the reduced density matrix of a subsystem,^., of, say, m qubits is defined as 
an average over the states of the other qubits: p\a{^ j) — Yli—i p(ik,jk). 

2.2 Quantum circuits with mixed states 

We describe the model of quantum circuits 0, [HJ, with mixed states Q: A quantum unitary gate of order k 
is a complex unitary matrix of size 2 fc x 2 fc . A density matrix p will transform by the gate to g o p = UpU\ 
where U is the extension of U . Using density matrices one can also define a non unitary gate: A measurement 
gate of order k is a complex hermitian matrix of size 2 k x 2 fe . A density matrix p will transform by the gate to 
g o p = O o (p). A Quantum circuit is a directed acyclic graph with n inputs and n outputs. Each node v in the 
graph is labeled by a quantum gate g v . The in-degree and out-degree of v are equal to the order of g v . Some of 
the outputs are labeled "result" to indicate that these are the qubits that will give the output of the circuit. The 
wires in the circuit correspond to qubits. An initial density matrix p transforms by a circuit Q to a final density 
matrix Q o p = g t o ... o g 2 o g l o p, where the gates <?t-..<?i are applied in a topological order. For an input string i, 
the initial density matrix is pu\. The output of the circuit is the outcome of applying basic measurements of the 
result qubits, on the final density matrix Q o pu\. Since the outcomes of measurements are random, the function 
that the circuit computes is a probabilistic function, i.e. for input i it outputs strings according to a distribution 
which depends on i. 

3 Quantum mediums with faults 

In this section we define quantum mediums subjected to probabilistic single qubit faults of one type. We then 
show that the probabilistic fault model is equivalent to a deterministic weak fault model and to a model with many 
possible types of faults. 

3.1 Quantum mediums with probabilistic faults 

In order to incorporate faults into quantum circuits, we need to add a clock, since the process of errors depends 
on time. A quantum circuit will be associated a timing, which indicates the time step each gate is applied. The 
timing should be consistent with the circuit: if a gate gets as input output of another gate, then it should be applied 
at least one time step after the other gate. In a faultless computation, one can assume without loss of generality 
that all the qubits are initialized and output together. In the presence of faults, it is still not clear whether allowing 
to input constants at any time helps the computation, since using qubits for a limited time might protect them 
against faults. Taking the stronger model, we permit also that qubits are input and output at different times, and 
we say a qubit is alive from t\ to t% if it is input to the circuit at t\ and output at ti. We will assume all timings 
start at 0, and denote by T the last time step. 

We would like to incorporate in the model single qubit faults. With out loss of generality, a single qubit fault, F, 
operating on the q'th qubit in the density matrix p of n qubits, can be described by adding a qubit to the system, 
that represents the environment, and letting a quantum gate to operate on it and the q'th qubit, and then taking 
the reduced density matrix to the original qubits. The fault is thus specified by the gate. The resulting density 
matrix is denoted as F q o p. In this model, a live qubit undergoes the fault F with some probability r\ per time 
step, and rj is called the decoherence rate. (One can allow the fault to depend on the time and on the exact qubit. 
All the results of this paper still hold.) 
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We now define the medium that realizes a quantum circuit with a certain timing, and is subjected to a certain 
dccoherence rate: 

definition 1 A quantum medium M(Q, t, 77}) is a triplet: A quantum circuit Q, a timing r on Q, and a decoherence 
rate < 77 < 1. 

We define a fault path, which is a possible list of where and when the faults occur: 

definition 2 Let M{Q,r,n) be a medium of n qubits, which lives T steps. Denote by t\,t\ the times of berth and 
death of the q'th qubit. A fault path is a set of different pairs (qj,tj), where 1 < qj < n, t\ 3 < tj < if ■ 

We associate with each path a weight, which is the probability for it to occur: 

definition 3 Let M{Q,t,w) be a medium of n qubits which has V pairs (q,t) where faults can occur. The weight 
of a path cr is w(a) = efH (1 - d) v ~^ . 

M computes a probabilistic function, where the distribution f(i) which M computes for an input i is defined 
as a weighted average over the outputs which the final density matrix of each fault path generates. To formalize 
this, we need notions of development in time of the medium. For a medium of n qubits, an evolution is a possible 
such development, i.e. a list {E(i)}f =0 of density matrices of n qubits. There are special evolutions which develop 
subjected to a certain fault path, and besides that "obey" the gates: 

definition 4 let M be a medium of n qubits, F a fault, a a fault path. An evolution E is called a "cr, F -trajectory 
of M" if it satisfies E(t + 1) = F at o FJ fe g\ o E(t) where g\ are the gates applied at time t, and F at are faults 
applied in a at time t. 

We can now define the function computed by the medium when subjected to faults F: 

definition 5 Let the medium M be subjected to faults F. For an input string i let pi = ^ CT w(a)Ef(T), where E? 
is an F, a -trajectory of M which satisfies EJ(0) = . The probability for M to output a string j for the input i 
is the probability to get an outcome j when applying basic measurements on the result qubits in pt. 

3.2 Equivalence to a weak interactions 

In the weak interactions model, each qubit undergoes a weak fault, deterministically, each time step. The weak 
fault operates on a density matrix in the following way: WF q o p = (1 — rj)p + r]F q o p. where F is a single qubit 
fault, and 77 is the strength of the interaction with the environment. The first term serves as a guard, so that the 
change in the density matrix is limited. To show that this model is equivalent to the probabilistic model, it suffices 
to show that there is no difference between a mixed state which is a weighted sum of mixed states, and a mixed 
state prepared by choosing, with the correct probability, one of the mixed states in the sum. 

lemma 1 Let p = p\p\ + P2P2- The mixed state p is quantumly equivalent to a mixed state prepared by choosing 
with probability pi the density matrix pi. 

Proof: The distinguish ability of quantum states is by measurements. By definition, the statistics of any mea- 
surement will be the same in the two cases. I 

3.3 Equivalence to a model with many possible faults 

We can also define a model which seems more physical: At each time step, and for each live qubit, different 
possible faults F i; can occur with probabilities 7^, where 77 = J2i Vi- K might seem surprising at first sight, but this 
model is equivalent to our probabilistic model since we can define the one fault in our model to be a probabilistic 
combination of all the faults: F o p — A 77^ o p. 
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4 Polynomial simulations of sequential mediums 

We describe a simulation of a quantum compuation by a quantum medium, on a randomized classical Turing 
machine. As described in the introduction, it suffices for our purposes to regard only collapse faults. We show 
that the expected slowdown of the simulation, for sequential mediums and mediums which operate gates only on 
0(log(n)) qubits each time step, is polynomial in the time and number of qubits of the medium, for any non-zero 
decoherence rate, rj > 0. 

To study the effects of decoherence on QTM one must refer to some possible physical realizations of QTMs. 
Our results apply immediately to the possible realization of QTM by sequential Quantum circuits, suggested by 
Yao[fL9||. At the end of this section we informally describe how to extend the results to general multi-head QTM. 

4.1 A classical simulation of quantum mediums 

We define an algorithm A(M, F)(i) which gets as an input a description of a Quantum medium M, a single 
qubit observable F, and an input string i, and outputs the function that the medium outputs when subjected to 
collapses of type F. A presents a density matrix, p, of n qubits, by a configuration: a list of density matrices, pj, 
each pj describing a set of entangled qubits, or a "cluster" , such that there is no entanglement between different 
clusters, p is the tensor product of the matrices in the configuration, let us describe the probabilistic algorithm A: 
Input: 

Description of a medium M = (Q, r, d), 
string of n qubits, i. 

The eigenstates |ei), |e2), of the observable F. 
Algorithm: 

- A chooses a decoherence path a with probability io(cr), according to the decoherence rate n. 

- A initializes the configurations conf(0),conf*(Q) to be a list of matrices of one qubit, |0)(0| or |1)(1| that 
correspond to i. 

- For t=l to T, A does: 

a) To simulate the computation step: Let g be a gate that operates at time t. Let p q i be the matrices in conf(t) 
that describe qubits on which g operates. A replaces them in the configuration by the matrix go (p qi <g) ...p qi ). For 
two qubits in the same cluster described by p q , we take p q in the tensor product only once. This is done for all the 
gates of the t'th time step to give conf*(t). 

b) To simulate a decoherence step: Let (g, t) be a pair in a. Let p q be the matrix in conf* (t) that describes the 
qubit q, and maybe a set B of more qubits. A computes, from p q , the probability for q to collapse on each of the 
eigenvalues of F. Than A tosses a coin with this bias, to decide on which of the eigenstates the qubit collapses. 
Let the result be |ei). A replaces p q in conf*(t) by the matrix |ei)(ei|, and the reduced matrix F ei o (p q )\B This 
simulates a measurement of the q'th qubit, conditioned that the outcome was e\. 

Output: To output the results of the computation, A imitates basic measurements of the result qubits, as in the 
decoherence step, and outputs the outcome. 

To show that A computes the same function as the quantum medium, we observe that the simulation actually 
computes the F, a — trajectory of the medium, for the path it chose, except one thing: The simulation makes a 
conditioned measurement, and continues according to the result, while the effect of the measurement faults is of 
an unconditioned measurement. However: 

lemma 2 The output distribution of A(M, F)(i) equals the distribution /m(«) that M computes when subjected to 
faults F , for an input i. 

Proof: By lemma |l|, the mixed state F q o p = Pi|ei)(ei| <g> F ei o p\ B + P2\e2) (e2 1 <8> F e2 o p\ B can not be distinguished 
from a mixed state prepared from | } (e^ | ® F ei o p\ B with probability pi. Therefore the effect of the unconditioned 
measurement is the same as taking with the correct probabilities the conditioned measurements. | 
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4.2 The simulation is efficient. 



To analyze the cost of the simulation, we divide the qubits each time step to two sets: The individual qubits, 
that in the configuration of A at time t, are described by their own density matrix, and the other qubits, which are 
entangled between themselves. Each individual qubit is described in the configuration by a density matrix with 4 
entries, but the description of the entangled qubits might be exponential in their number. Denote by K{t), K*(t) 
the sizes of sets of the non-individual qubits at time t in conf(t), conf*{t) in the simulation. To bound the cost of 
the simulation from above, we first prove that the probability for the number of non-individual qubits to be large 
is exponentially small. 

lemma 3 Let M be a medium which operates on no more than P qubits each step, subjected to collapses of type 
F, with n > 0. Consider A(M, F)(i) . There exist constants a > 4, and b > such that Pr(K(t) > bP + j) < -p-, 
Pr(K*(t)>bP + j + P)<±. 

Proof: We use induction on t. For t = 0, K(0) — 0, K*(0) < P. For the induction step, assume Pr(K(t) > 
bP + j) < ~ smdPr(K*(t) > bP + j + P) < i fort. To analyze Pr(K(t + 1)), We apply decoherence on K*(t). 
The probability to reduce the number of non-individual qubits is binomial so: 

oo 

Pr(K(t + l)>bP + i) = Y J Pr(K(t + 1) = bP + i + j) = 



3,1=0 



3=0 




a — 1 \ a — r\ ) \a — rj / 

We want this to be smaller than ~, for any i. This implies two demands: (1) < ^ which is achieved for 
a > 1. (2) The coefficient of the geometrical sequence must be smaller then 1. There exists b that fixes that since 
^a-i" 1S strictly smaller than 1, because n > 0. 

For Pr(K*(t + l)) we have that K*(t + 1) < K(t+l) + P. Then Pr(K*(t+l) >bP + i + P)< Pr(K(t + l) > 
bP + i) <i.| 

We can now prove that the expected number of entries of density matrices that A writes is polynomial: 

lemma 4 Let M(Q,T,r[) be a medium of n qubits with parallelism 0{log(n)), n > 0, F a measurement fault, i an 
input string. Consider A(M, F)(i). The expected number of entries of the density matrices in all the configurations 
that A writes is polynomial in n and T. 

Proof: We use the definitions of the constants a and b from lemma|. We compute the expectation value of m(t), 
the number of entries for the matrices of the K(t) non-individual qubits, at time t: E(m(t)) = J2m=i P r ( m (t) — 
m)m < Yl'jLo P r (K(t) = < Y^j=o 1 4 J + J2hLo i r ^ bP+ '' wnere the last inequality is by lemma |[ The sum 
is finite since a > 4. The expectation is constant if the parallelism, P, is constant, and for P = clog(n), the 
expectation is 8(n 2cb ). At time t we also have to write 4 entries for each small density matrix, which adds at most 
An entries. The expected total number of entries, (during the whole simulation) is the sum of the expectations in 
each time step. Thus it is polynomial in n and linear in T . | 

Theorem 1 Let M be a medium of n qubits, life length T , with n > 0, and 0{log{n)) parallelism, F a measurement 
fault. The expected cost of A(M, F)(i) will be polynomial in T and n. 

Proof: By lemmas 0,0 and the fact that one entry of a density matrix can be written in poly(n,T) bits, if the 
specification of the gates is by poly(n,T) bits. | 
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4.3 Simulating quantum Turing machines 



We briefly mention how to extend the above result to other possible models of QTM. We assume a decoher- 
ence model for multi-head QTM where besides single site faults, we allow faults on the heads' positions and the 
processor's state. Again we refer to the model of random collapses, where we can show: 

Theorem 2 Let Q be a QTM, subjected to measurement faults F of the cites, the heads' position and the machine's 
state, with r\ > 0. A classical probabilistic Turing machine can simulate T computation steps of Q on inputs of 
size n in time polynomial in T and n. 

Proof: Sketch. After every constant number of steps, the head's position will collapse to one place. Since at 
each time step it can only triple the number of sites it works on (by going Left, right or staying in place) the 
expected number of different paths of the head at each moment will remain constant. The simulation will remain 
efficient when the number of heads is O(logn) since then the number of paths is polynomial. 



5 Simulating mediums with general parallelism 



In this section we study the effects of decoherence on the computational power of general Quantum mediums, 
by applying the simulation A, defined in section 4.1, for mediums with general parallelism. We assume gates of 
fan- in < 2, but the results hold for any constant fan-in. We show analytically that the cost of the simulation has 
an expected polynomial slowdown for any medium with high enough decoherence. If the medium uses gates on 
more then one qubit not too frequently, the simulation is still polynomial even for lower decoherence rates. The 
simulation is proved to be exponential for a general random medium with low enough decoherence rate. Finally we 
present computer experiments, showing a phase transition in the expected cost at some decoherence rate, where 
the cost "jumps" from polynomial to exponential. These transition points are found for random circuits and 
nearest neighbor circuits. Known results ]T(| imply that for random circuits the transition happens in a range of 
decoherence rates which is 0( / og („) )■ For any sequence of quantum circuits for which there is a gap between the 
costs for high and low decohcrcnces, the transition happens in a range which goes to zero as the size of the circuit 
goes to infinity. 



5.1 Mediums with high decoherence rate. 



It is proved that for high enough decoherence, the expected slowdown of the simulation A is polynomial. We 
need a more refined analysis than in the sequential case, since the line of reasoning applied there leads in the 
general parallelism case to a trivial exponential upper bound on the cost. Recall that a cluster in a configuration 
of A is a set of qubits described by the same density matrix in the configuration. In the analysis of the sequential 
case we referred to the total number of non-individual qubits, i.e qubits which are in clusters of size bigger than 
one. Here we take a closer look and consider the sizes of each cluster alone. We define K q (t),K*(t) as the size of 
q's cluster in conf(t),conf*(t), respectively. As long as the sizes of the clusters are 0(log(n)), the simulation is 
polynomial in n, since the number of entries in each density matrix in conf(t) will be polynomial. 

We intend to show by induction on t that the probability for K q (t) to be large decays exponentially. The 
difficulty here is to bound the change in the clusters' sizes during the computation step. If we could bound this, 
so that the clusters do not grow too much during the computation step, and if the decoherence is high enough 
to shrink them back down, then the induction step will work. In the following we analyze separately the change 
in the distribution of sizes of clusters, during the computation step, and during the decoherence step. Then we 
combine the two to show that the cost of the simulation is polynomial for high enough decoherence. 



5.1.1 Clusters' Growth in a computation step 

Assuming an exponential decaying bound on the sizes of clusters before the computation step, we can give an 
exponential decaying bound on the sizes of clusters after the computation step. 

To do that we observe a connection to a branching process. Let the clusters in the configuration at time t, 
conf(t), be A, B, C, ... Let M be the next computational step. The question is what is the cluster's size of a qubit 
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a <E A after M is applied. Recall that we assume gates of fan-in two at the most. If a gate in M operates on two 
qubits then they are called mates. We associate with a a tree of which a is the root, a's sons in the tree are the 
qubits in A. Suppose a\ is a son of a and a mate of b S B. Then all members in B are the sons of ai, and so on. 
In each stage we ignore qubits that already appear in the tree. The distribution on sizes of clusters corresponds 
to a distribution on the number of sons in a branching process. The distribution on the sizes of clusters after the 
computation corresponds to the sizes of trees in the branching process. If a has no mate, K*(t), which is the size 
of a's cluster after M is applied, equals the size of a's tree. If a has a mate, c, then K*(t) is the sum of sizes of a's 
tree and of c's tree, where again duplicated qubits are ignored. 

We first consider a simplified version of the above process where (1) All gates in M are of fan-in 2, and (2) Gates 
in M operate only on qubits from different clusters. With these simplifications, the process that generates a tree, 
defined above, turns out to be a simple branching process, for which we can compute the probability distribution 
over the total size of two trees, given the distribution on the number of sons. 

lemma 5 Let B be a branching process, with Pr(z = j) < ^j,a > 8 for z, the number of sons. Let T a ,Tt,, be the 
sizes of two independent trees generated by B. Let L — 2T a + 2Tb — 2. Then Pr(L — i) < a(y/2 — ^-) 2 (^J^T- 

Proof: For the standard proofs in the following discussion, consult |fl3| . Let the generating functions^ of z, T, L 
be G z , Gt, Gl- In the following formulas, the first jl3) connects the generating functions of the size of a tree T to 
the number of sons Z By definition of L we have. The second we have by definition of L: 

G L+2 (x) = (G T (x 2 )) 2 , Gr(x) = xG z (G T (x)) (1) 

We don't know the generating function for the variable z since we only have a bound on the distribution. Let 
us refer to the exponential decaying bound as a pseudo distribution, although it is not normalized to one, and 
define "pseudo generating functions", (denoted by a tilde), such that G z {x) will "generate" the exponential decay 
i: G z (x) = TZoW xi = 5=5> and ^t(x), G l (x) will satisfy equations |: G T (x) = xG z (G T {x)) , G L+2 {x) = 
(Gt(x 2 )) 2 . The coefficients Zi,ti,k in the analytic expansion of Gz{x) Gt{x), and Gl+2(x) are smaller then the 
corresponding coefficients Zi,ti,U of the "pseudo generating function," when the expansions exist. To prove this, 
note that the assumptions of the lemma imply Zi < 2j. Writing equation [l] using analytic expansions, one can show 
by induction that t\ < t\. Using this and writing equation |l| in analytic expansion, we have h < h- 

We intend to give an upper bound on Zj, since this will imply an upper bound on li + 2 = Pr(L = i) and therefore 
on Pr(L > i). Given G z (x), we can solve Gt{x) = xG z (Gt(x)) for Gt{x), which gives a quadratic equation. 
We choose the "minus" solution of the quadratic equation, which corresponds to the correct choice in the case of 
G T {x). It gives G T (f ) = (x/ ^ 1)a , which implies 

iV ~ 2 - 1)2a2 =(G T £)) 2 = G L+2 (S- (2) 



Gt(x) is analytic for |x| < |, so Gl +2 (j/) is analytic at y = Vg, and we can write: { ~% = Si*lo^(\/T) 
which gives k < (\/2 — 1 ) 2 fyi) implies the desired result. | 



We use the connection to the branching process defined above, to show an upper bound on the clusters' size 
after the computation step, given the upper bound before the computation step. 



lemma 6 Let a > 8, and C 



(V2-lf 7 7 =. Pr(k q (t) >j)<^r implies Pr(K*(t) > i) < C(^Y 



Proof: Compare our process to the simplified branching process, where identifying k q {t) — 1 with the variable z 
and k*(t) with L. Our process differs in two aspects: It might connect mates inside a cluster, and a qubit might 
not have a mate at all, i.e might not have a chance to give birth to children. The only effect of these fact can be 
to reduce the final cluster. | 

4 The generating function Ga(%) of a random positive integer A, is the expectation of x A : E(x A ) = J^°l Pr(A = i)x* 
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5.2 The decoherence step 



We now show that with strong enough decoherence rate, the distribution over sizes of clusters is "pulled" by 
decoherence back below the original exponential decaying bound: 



lemma 7 3 <q x < 1, a > 8 such that V rj > rji Pr(K*(t) > i) < C(^/f) 1 with C = (VS - l) 2 yf^/g , impHes 
Pr(fe(t+l)>0<3^r. 

Proof: We find 771 < 0.97. We have to calculate how a distribution over sizes of clusters changes because of the 
decoherence. The probability to reduce the number of qubits in a cluster is binomial so: 

00 

Pr(k(t + 1) > I) = Pr(k(t + l) = l + r ) = 

r=0 

00 00 / 7 \ 

££Pr(**(*)=i + r+p)( l + r + P ) V P(l- V ) l+r < 



oo 00 / . \ 

E^(^(*)>«+p)E( P +p )v p (i-v) l+r < 

F~Q~ 00 00 Ao" / 7 

r— p * 



(1 - 



r=0 (1-^I)'+H-1 

We want that for all I Pr(k(t + 1) > /) < a(-) , so we demand that the factor of the geometrical sequence is 
smaller than - and that the first term is smaller then a. The demands are satisfied with a = 30 and ij > 0.97.1 

5.2.1 Closing the proof 

Theorem 3 There exists a constant rji < 1 such that for any medium M on n qubits, with decoherence rate 
rj > T)i, which lives for T steps, and any measurement fault F and input string i, the expected cost of A{M, F)(i) 
is polynomial in n and T . 

proof: We show that for A of mediums with a decoherence rate higher than 771 defined in lemma M, and 

for any qubit q, the sizes of q's clusters in the algorithm A satisfy for some a > 8 Pr(k q (t) > j) < -j and 
l — 3 

Pr(K*(t) >i)< C(^/f Y for C = {s/2- l) 2 ^ 2 ^ . This is done by induction on t. The basis of the induction is 

trivial, as K q (0) = l,i^*(0) = 1 always. The induction step is done in two stages, by lemmas |6| and [?|. With this 
exponential decaying bound on the distribution over sizes of clusters, we can proceed to show that the expected 
number of entries of matrices that A writes is polynomial, as in lemma ^, and since each entry can be written 
using a polynomial number of bits, ( if the gates are specified in polynomial number,) this completes the proof, 
with 771 = 0.97.1 

The numerical lower bound of r/i — 0.97 is ofcource not a reasonable physical number, to even start computation 
with; It is high because it is good for any quantum circuit. Most circuits are simulated efficiently with much lower 



decoherence rate(see section |5.5|) 
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5.3 The importance of the frequency of many-qubits-gates 

A quantum circuit that uses single qubit gates often, is simulated efficiently already for low decoherence rates, 
since single qubit gates do not participate in the struggle against decoherence. To see that, consider a quantum 
circuit that applies a many-qubits gate on each qubit once every v time steps at the most. One can treat these v 
time steps as a one time step, in our branching process analysis. During this "extended" time step, the decoherence 
had it's chance to operate v times. The effective decoherence rate, i.e. the rate of decoherence with respect to the 
extended time step, is therefore 1 — (1 — rf) v ^> 77. Therefore already for very small decoherence rates rj the effective 
decoherence will be high enough for the simulation to be efficient. 

5.4 The simulation is exponential for low decoherence in random circuits 

In this subsection we turn to the other side of the decoherence rate scale, and ask what is the expected cost for 
low decoherence rates. Off cource, there are circuits which can be simulated efficiently even with no decoherence 
at all. However, for random circuits the expected cost of the simulation is exponential for 77 < 0.5. By a random 
circuit we mean a circuit in which each time step a random matching between qubits is chosen, and this matching 
corresponds to gates of fan- in two operating on the qubits. 

lemma 8 The simulation of a random circuit subjected to decoherence rate n < 0.5 is exponential. 

Proof: We prove that when the decoherence rate is 77 < 0.5 there is a cluster of size linear in n. During one 
computation step a cluster of an qubits, will grow to (a + a(l — a))n since the probability for qubits in the 
cluster to be matched outside the cluster is (1 — a), and these matching add at least one qubit to the cluster. 
The decoherence step multiplies this by 77. Hence the original cluster size is multiplied each time step by a factor 
(2 — a)r). For 77 > 0.5 this is bigger than 1 for small a, and the cluster will grow until it is of size (2 — which 

is linear in n. The rate of growth of the cluster is exponential, because as long as a is smaller than 2 — the 
factor is bigger than 1 + e. 

5.5 Phase transitions: experimental results 

We have performed computer experiments to investigate how the clusters sizes depend on the decoherence rate. 
To study the random algorithm case, we started with n clusters of size one, and each time step chose a random 
matching between qubits. A match represented a gate on two qubits, so their clusters were joined. To imitate 
faults, we randomly and independently separated a qubit from it's cluster with probability 77, each time step. The 
system exhibits a phase transition at decoherence 770 ~ 0.64. For 77 > 770 the clusters' sizes stay smaller than log(n), 
but for 77 < 770 a "giant" cluster of size linear in n appears, which means that the density matrix describing it is 
exponential. 

We also studied experimentally a one dimensional array of qubits, where gates where allowed to operate on 
nearest neighbors only, and each qubit was matched with it's left neighbor and right neighbor alternately. We 
found the transition point at 77 ss 0.5. 

6 Conclusions 

We feel that this work gives insight on the effects of decoherence on quantum computation, and emphasizes the 
importance of 0{n) parallelism in a quantum computation subjected to decoherence processes. 

Our results with respect to sequential quantum computers, and quantum computers with (log(n)) parallelism 
show that any amount of decoherence destroys all the quantum computational power. 

Considering parallel computers, an upper bound was given on the decoherence rate which any quantum computer 
can hope to tolerate. An important role is given to the size of the fan-in used in the computation, and to the 
frequency of many-qubits-gates. Computers with many single qubit gates seem more sensitive to decoherence, 
and can be simulated efficiently even with very low decoherence rates. This should be considered when planning 
implementations of quantum gates by low- fan-in gates H |l). It might also be interesting to quantify, according to 
these observations, the stability of different algorithms against decoherence. 
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DiVincenzo |9| raises the question of whether the time evolution of an open quantum system is a powerful 
computational tool. Our conclusions are that at least within the model we have used, the answer is negative. 

An open question is whether a stronger simulation than ours exists, that is efficient for lower decoherence rates 
or even with no decoherence at all. 

Early results Jl5|, |j| indicate the possibility of an efficient simulation of an ideal QC by quantum mediums 
subjected to low decoherence rates. Though our results might seem as a step towards negating the fascinating 
possibility of a physical realization of quantum computers, we actually view these results as optimistic, since phase 
transitions are known to be robust, which implies that maybe quantum computers with low decoherence rates 
indeed can not be simulated efficiently. 

7 Acknowledgments 

We wish to thank Noam Nisan, Nati Linial, Yuval Peres, Ehud Fricdgut, Amnon Ta-Shma, David DiVincenzo, 
Arthur Ekert, Avner Magen for very helpful conversations, essential remarks and corrections. A special thank to 
Yoram Cohen for invaluable help in so many aspects. This work was initiated during the Quantum Computation 
Workshop at Torino, supported by the ISI Foundation. 

References 

[1] A.Barenco, C.H. Bennett, R. Cleve, D. P. DiVincenzo, N. Margolus, P. Shor, T. Sleator, J. Smolin, and H.Weinfurter. 
Elementary gates for quantum computation. In Phys.Rev.A 52, pages 3457-3467, 1995. 

[2] D. Aharonov and N. Nisan. Quantum circuits with mixed states, in preparation. 

[3] A.R.Calderbank and P.W.Shor. Good quantum error correcting codes exist. In phys.Rev.A, to appear, pages quant- 
ph/9512032, 1995. 

[4] E. Bernstein and U. Vazirani. Quantum complexity theory. In Proceedings of the Twenty-Fifth Annual ACM Symposium 
on the Theory of Computing, pages 11-20, 1993. 

[5] I. L. Chuang, R. Laflamme, P. W. Shor, and W. H. Zurek. Quantum computers, factoring and decoherence. Science, 
270:1633-1635, 1995. 

[6] D. DiVincenzo. Two-bit gates are universal for quantum computation. 1995. 

[7] D. Deutch. Quantum networks. In Proc. Roy. Soc. Lond, Vol. A400, 1989. 

[8] D. Deutsch. Quantum theory, the church-turing principle and the universal quantum computer. In Proc. Roy. Soc. 
Lond, Vol. A400, pages 96-117, 1985. 

[9] DiVincenzo. Quantum computation. Science, 270, 1995. 

[10] E.Friedgut and G. Kalai. Every monotone graph property has a sharp threshold, to appear. 

[11] J.J.Saqurai. Modern Quantum Mechanics, revised edition. Addison Wesley, 1994. 

[12] S. Lloyd. A potentially realizable computer. In Science, Vol 261, pages 1569-1571, 1993. 

[13] N.Alon and Spencer. The probabilistic method. 

[14] G. M. Palma, K.A. Suominen, and A. Ekert. Quantum computation and dissipation. Proc. Roy. Soc. Lond., 1995. 

[15] P. W. Shor. Fault-tolerant quantum computation. In 37th Annual Symposium on Foundations of Computer Science, 
1996. 

[16] P.W. Shor. Algorithms for quantum computation: Discrete logarithms and factoring. In 35th Annual Symposium on 
Foundations of Computer Science, pages 124-134, 1994. 

[17] W. G. Unruh. Maintaining coherence in quantum computers. Technical report, University of Vancouver, 1994. quant- 
ph/9406058. 

[18] W.H. Zurek. Decoherence and the transition from quantum to classical. In Physics today 44(10), pages 36-44, 1991. 

[19] A. Yao. Quantum circuit complexity. In 34th Annual Symposium on Foundations of Computer Science, pages 352-361, 
1993. 



12 



